The classical skeleton of open quantum chaotic maps 
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Abstract 

We have studied two complementary decoherence measures purity and fidelity for a generic diffusive noise in two different 
chaotic systems (the baker and the cat maps). For both quantities, we have found classical structures in quantum mechanics 
- the scar functions - that are specially stable when subjected to environmental perturbations. We show that these quantum 
states constructed on classical invariants are the most robust significant quantum distributions in generic dissipative maps. 



PACS numbers: 05.45.Mt, 03.65.Sq, 03.65.Yz 



I. INTRODUCTION 

In quantum chaos (i.e. the study of quantum systems 
whose classical counterparts are chaotic), the complete 
description of the eigenfunctions is still an open prob- 
lem. While there is a semiclassical method for obtaining 
the eigenstates of an integrable system (the well-known 
EBK/WKB quantization scheme [1, 2]), the same prob- 
lem in the chaotic case has proven harder to solve. One 
of the major advances in this sense has been Gutzwiller's 
theory of periodic orbits for the quantum chaotic spectra 
[3], devised in the early 70's. 

A recently developed approach allows to obtain all the 
information of a generic quantum system by just using 
the shortest periodic orbits of its classical counterpart 
[4]. This formalism has been successfully applied to the 
description of the eigenstates of many chaotic quantum 
systems like for example, the Bunimovich's billiard [5], 
the cat map [6] and the baker map [7] . These results sug- 
gest that states built using only this classical information 
(the so-called scar functions) constitute the skeleton of 
the eigenstates of any quantum chaotic system. 

On the other hand, taking into account the restora- 
tion of the classical dynamics induced by decoherence in 
an open quantum system [8] and the behavior of distri- 
butions against external perturbations [9- ] we can ask 
ourselves: are scar functions, as a consequence of their 
classical content, more robust than other states facing an 
environment- induced decoherence? Trying to answer this 
question, in [12] we have studied a very simple model of 
a chaotic quantum system interacting with an environ- 
ment that produces dissipation and decoherence (an open 
quantum map). In this article, we exhaustively expand 
the results presented there by studying different open 
maps. This leads us to our main result, i.e. we show 
here that scar functions are the most robust significant 
quantum distributions in generic dissipative maps. 

In this paper we study the behavior of classically moti- 
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vated states (by analyzing two different measures, namely 
the purity and fidelity) corresponding to two paradig- 
matic systems in quantum chaos, the baker map and the 
cat map on the torus [13, 14]. We introduce decoherence 
by means of a diffusive noise model [15]. We have orga- 
nized this work as follows. In section 2, we describe the 
main theoretical tools we need throughout our investiga- 
tions (open quantum maps and scar functions). In Sec- 
tion 3 we define the maps considered. Section 4 presents 
our results. In section 5 we come to some conclusions. 



II. THEORETICAL TOOLS 

Open quantum maps are the simplest systems that 
capture all the essential features of chaotic dynamics 
and dissipation. As such, they are the ideal testbed for 
studying the effect of the environment on quantum dis- 
tributions. Our main interest is focused on determining 
how the classical information emerges from the quantum 
structures. For that purpose we will make use use of the 
scar functions, whose construction we explain at the end 
of this Section. 



A. Open quantum maps 

The quantization of maps on a compact phase space 
proceeds in two stages: a kinematic one, which estab- 
lishes the Hilbert space appropriate to the phase space 
geometry, and a dynamical one, which consists in defining 
a suitable quantum operator corresponding to the classi- 
cal dynamics. Finally, Kraus operators take into account 
the effect of the environment. 



1. Kinematics 

In this work, the classical phase space associated with 
the systems under investigation is the 2-dimensional 
torus T 2 = R 2 /Z 2 , consisting of a square of unit side 
with opposite sides identified. Points in this space have 
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coordinates (q,p) G [0, 1) x [0, 1). To be compatible with 
the phase space geometry, quantum wave functions must 
be periodic in both position and momentum (up to a 
phase) 



ip(q + l) = exp(-i27rx (? ) ip(q) 
= exp(i27TX P ) $(p) 



with 



dq e~^ qp i/j(q). 



(1) 
(2) 



(3) 



The phases 2i\\q^Xp are called Floquet's angles, with 
< Xqi Xp < 1- Periodicity in q and p implies 



2nhN 



1 



(4) 



This means that the Hilbert space of wave functions 1~Ln, 
the quantum counterpart of the classical compact phase 
space, is effectively finite-dimensional with dimension TV. 
In this context, the semi-classical limit h — >• is equiva- 
lent to taking TV — >• oo. 

The usual canonical commutation relations between 
position and momentum operators (which are the gen- 
erators of infinitesimal translations in phase space) do 
not hold in the finite-dimensional case, so we can't de- 
fine position and momentum operators q,p. However, we 
can define finite displacement operators U and V, whose 
form is analogous to that of the infinite-dimensional case. 
These operators are unitary and its eigenvectors form a 
basis for Hat [16-19]. 

The basis of position vectors \qj) for this space will be 
defined from the eigenvectors of the momentum displace- 
ment operator V 



V\qj) = exp 



27ri 



U + Xq) 



fa), je[0,N-l] (5) 



The position in phase space associated with \qj) is 

3 +Xq 



TV ' 



(6) 



Analogously, the momentum vectors \pk) satisfy 



U\p k ) = exp 



27ri 



(k + x P ) 



\pk) 



Xp 



N 



, ke[o,N-i\ 



(8) 



where U is the position displacement operator. 

These bases are related by means of a discrete Fourier 
transform 



(pk\qj) = N 2 exp 



■jlJ +X<z)(fc + Xp) 



= F?> k 



N 



(9) 



From the previous equations, it can be shown that [16- 
19] 

Ufa) = \q j+ i) (10) 

U N \ qj ) = \q j+N )=exp(-i27r Xq )\q j } (11) 

V\p k ) = \p k+1 ) (12) 

V N \pk) = \Pk+N) = exp (i27TX P ) \pk) (13) 
These operators satisfy the relation 



(14) 



With them, we can define a discrete version of the phase 
space displacement operator 



4, = -^expfe/c 



(15) 



with the property T] k = T_j _ fc . The set \ T jik \ of 

finite displacements forms a basis for the Hilbert space 
7-L N 2 = % N (g) 1-L* N of linear operators on Hn (Liouville 
space) with the Hilbert-Schmidt inner product 

=Tr(i t s) (16) 

This basis of N 2 displacement operators satisfies 

Tr(7;;.,7V./,-) <\,..r^Jy (17) 
so it constitutes a complete orthonormal set. 

2. Dynamics 

By virtue of the finite dimension of 'Hat, the quantum 
dynamics is given by a unitary TV x TV matrix Un and the 
system state evolution is obtained by means of a straight- 
forward matrix multiplication. In our case, this matrix 
will be the quantization of a classical (chaotic) map P on 
T 2 . This means that, given the mapping P(q,p), there 
exists a sequence of unitary operators Un (called quan- 
tum map) acting on Hn such that the so-called Egorov 
property is fulfilled, i.e. 



lim 

iV-^oo 



U^Op(f)U N -Op(foP) =0 V/eC7°°(T 2 ) 

(18) 

Op(/) represents the Weyl quantization of the observable 
/ [16]. 

3. Environment 

If the quantum system interacts with an environment, 
the elements (kets) of Hat no longer represent its state. 
In this situation, all that can be said about the system at 
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time t is encoded in its associated density operator p t [20- 
22]. The evolution of p t is given, under sufficiently gen- 
eral conditions, by a completely positive, trace-preserving 
map S of density matrices into density matrices called 
superoperator or quantum operation [21-23]. 

In this context, we define an open quantum map S 
[15, 24] as a map whose action can be written in the 
form of a product of two superoperators 

Pt+ 1 = S(pt) = D e M(p t ) (19) 

M(p t ) = Mp t M^ is a map that generates the unitary 
evolution of the system, M being the evolution operator 
that acts upon elements of the Hilbert space associated to 
the non-interacting system. D £ is a superoperator that 
models the interaction between system and environment 
according to a set of parameters e related to the specific 
type of interaction. This last superoperator is responsible 
for introducing noise in the -otherwise unitary- system 
evolution. 

In the present work, M will be the quantization of a 
classically chaotic map acting on the torus. In particu- 
lar, we will consider two completely chaotic maps: the 
baker map and the cat map [13, 14]. Besides, the noise 
superoperator will be expressed according to the Kraus 
representation [21-23], which in general can be expressed 
as 

7V 2 -1 

D £ (p t ) = KiPtKj (20) 

i=0 

with Y^i-K\Ki = 1 in order to preserve the trace of p t . 
The way the system interacts with the environment will 
be completely determined by the Kraus operators K{. 
In the following Section, we will define these operators 
explicitly. An example of the action of the environment 
in phase space can be seen in Figs. 1 and 2. 

B. Scar functions 

The semiclassical theory of short periodic orbits [4] is a 
formalism that allows to obtain all the quantum informa- 
tion of a chaotic Hamiltonian system in terms of a very 
small number of short periodic orbits. The main elements 
in this theory are the so-called scar functions. These 
are wavefunctions highly localized in the neighborhood 
of the classical periodic orbits and on their stable and 
unstable manifolds, satisfying a Bohr-Sommerfeld quan- 
tization condition along the trajectory. They are defined 
for Hamiltonian flows as 




where H is the system's Hamiltonian, T is of the order of 
Ehrenfest's time, and \4>t u be) is a wavefunction localized 
on the periodic orbit with Bohr-quantized energy Ebs- 



In [6, 7] the formalism has been adapted to quantum 
maps on the torus, and the resulting formula for scar 
functions is given in terms of a sum 

icr ) = E cos (H ) eiEBSt u^poZ) (22) 

where U is the evolution operator of the quantum map 
and |^po^) (called Periodic Orbit Mode or POM) is 
a sum of coherent states on the torus centered at the 
fixed points of a given periodic orbit, each one having a 
phase. In this case, the Ehrenfest time is T = A is 
the Lyapunov exponent of the map and N is the Hilbert 
space dimension. 

As an example, the upper right panels of Figs. 3 
and 4 show Husimi representations of scar functions con- 
structed for the baker map and for the cat map, respec- 
tively. It's clearly visible the enhancement of probability 
that these wavefunctions have on the corresponding pe- 
riodic orbit and its stable and unstable manifolds. 



III. SYSTEMS 
A. Baker map 

The first model of chaotic dynamics we consider is the 
baker map B : T 2 T 2 , given by the transformation 

[13, 14] 

(q',p') = B(q : p) = [2q-l2q\,(p+ [2q\)/2] (23) 

where |_^J stands for the integer part of x. This trans- 
formation is an area-preserving, uniformly hyperbolic, 
piecewise-linear and invertible map with Lyapunov expo- 
nent A = In 2. The vertical (horizontal) lines q = qo (p = 
Po) represent the stable (unstable) manifolds. 

The phase space has a very simple Markov partition 
consisting of two regions (q < 1/2 and q > 1/2) as- 
sociated with the symbols and 1, for which there is 
a complete symbolic dynamics. The action of the map 
upon symbols can be understood by means of the binary 
expansion of the coordinates 

(p\q) = ...i/_i- v Q vi ... — >• (pf\qf) = ... v-ivo ■ v x . . . 

(24) 

where q = Y^bLo ^2 _ ^ +1 ^ and p = J2 i J*!_ 1 V{2 % ' . Then, a 
periodic orbit of period L can be represented by a binary 
string v of length L. The coordinates of the first trajec- 
tory point (qo,po) on the periodic orbit can be obtained 
explicitly in terms of the binary string as qo = vw . . . = 
vj(2 L - 1) and p = v^v^v^ ... = ^/(2 L - 1), where 
v is the integer value of the string v which represents a 
binary number, and is the string formed by all L bits 
of v in reverse order. The other trajectory points can 
be easily calculated by iterations of the map or by cyclic 
shifts of v. 
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The unitary operator M that performs the closed 
quantum evolution is given in position representation by 

[25, 26] 



M = F 



N 



N/2 



N/2 



(25) 



where Fn is the TV-dimensional Fourier transform opera- 
tor whose matrix elements were defined in (3). Through- 
out the paper we assume for the quantum baker map 
a phase space with anti-symmetric boundary conditions 
(Xq — Xp — 1/2) in order to preserve the classical map 
symmetries [26]. 



B. Cat map 

Another simple model with strongly chaotic dynamics 
on T 2 is the cat map [13, 14]. It's an invertible, area- 
preserving canonical transformation A whose matrix has 
integer entries, and with Tr(A) > 2 to ensure hyperbol- 
icity. A common choice for A is 



2 1 

3 2 



mod 1 



(26) 



The Lyapunov exponent for this map is A = In (2 + 
\/3) ~ 1.317. The expanding and contracting eigenspaces 
through the origin are given by £ u = (— >/3, 1) and 
£ s = (>/3, 1). The irrational slope of the two directions 
implies that stable and unstable linear manifolds are 
densely distributed over the torus. 

The map is quantized by means of its generating func- 
tion [14, 27], giving a unitary propagator M whose matrix 
elements in position representation are 



1 



Noise model 



exp 



2?ri , 9 



jk 



(27) 



We define the superoperator D £ by means of transla- 
tion operators on the torus [15] 



De(Pt) 



N-l 

V" Ce {j,k) f jtk p t f} 



j,k=0 



(28) 



. These translation operators (our Kraus operators) are 
defined by Eq. (15). To preserve Tr (pt) we assume 
^fk^o^U^) — 1 • This way, the coefficient c € (j,k) 
represents the probability of a translation being applied 
on the system in the direction (j, k). Defining this func- 
tion as a periodized Gaussian (to satisfy the boundary 
conditions on the torus) 



fi J u— — oo 



(j - iiN) 2 + (k- vNf 

H€) 2 



(29) 



we obtain a noise superoperator which has the effect of 
diffusing the state on a region of radius « e in phase 
space. The consequence of this incoherent superposition 
of translations is decoherence, which can be visualized 
as the suppression of the small scale interference fringes 
in the Wigner representation. After a short time, the 
Wigner function becomes positive and the state appears 
"smeared out" in phase space. As previously mentioned, 
this behavior can be seen in Figs. 1 and 2, which show the 
effect of noise over the discrete Wigner function [16, 17, 
19] of a superposition of coherent states. The parameter 
e can be interpreted as a measure of the coupling between 
the system and the environment. 



IV. RESULTS 

In order to quantify the stability of the states of inter- 
est against decoherence, we have calculated the purity 



P(t) = Tr(^) 
and the fidelity or autocorrelation function 



F(t) = VMptl 



(30) 



(31) 



as functions of time. In (31), 1^) represents the pure 
initial state, hence po = 

Purity is a measure of the correlation degree between 
the system and the environment, and its evolution in time 
indicates how fast the system loses coherence. Fidelity 
can be interpreted as the distance between the evolved 
state and the initial state. Its development in time al- 
lows us to measure the velocity with which the evolved 
state "moves away" from the initial one under the action 
of the noisy dynamics. Complement arily, this difference 
between the initial and evolved states can be observed in 
terms of a density operator representation in phase space, 
like Husimi or Wigner distribution functions [16, 28]. 

We have studied the evolution of these quantities from 
initial states defined as scar functions, periodic orbit 
modes and eigenstates of the unitary quantum map. For 
each system, we have built the scar functions and the 
periodic orbit modes on different short periodic orbits of 
the corresponding classical map (without noise) and then 
compared their evolution to those of the map eigenstates, 
in particular with those localized on or nearer the same 
orbits. We have also studied the effect of varying the cou- 
pling with the environment by means of taking different 
values for the parameter e. 

In the following Figures, we show some typical results 
of our numerical calculations. They are illustrative of an 
exhaustive exploration of the eigenstates of our systems. 
Fig. 3 shows the behavior of purity and fidelity for the 
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Figure 1: Action of the noise model in phase space. The first panel (left) shows the discrete Wigner function of a 
superposition of two coherent states centered at (0.35,0.35) and (0.65,0.65). The value of the Wigner function is 
shown using a gray scale from white (minimum, negative) to black (maximum, positive) (same scale for all panels). 
As time goes on (from left to right) the noise superoperator D £ acting on the state washes out the interference 

fringes. We have taken e = 0.05 and N = 100. 




Figure 2: Action of the noise superoperator on the same initial state as in Fig. 1, but with stronger coupling 
{e = 0.1). In this case the interference fringes disappear already in the first application of the noise superoperator. 
The gray scale chosen in this case enhances the visualization of the Wigner function spreading in phase space. 



case of the noisy baker map model. The states are lo- 
calized over a period 2 orbit with symbolic code 01, as 
can be seen in the Husimi representations of the upper 
panel. The overlap between the scar function and the 
map eigenstate is 0.828, so the states are quite similar 
in terms of this distance measure. However, there is a 
visibly different behavior in terms of purity and fidelity, 
favoring the semiclassically constructed states in general, 
and the scar functions in particular, which lose purity and 
fidelity at a slower pace. 

For the noisy cat map we have a similar behavior, but 
the differences between states in terms of purity and fi- 
delity evolution are less pronounced, as can be seen in 
Fig. 4. In this case, the scarring of map eigenstates by 
periodic orbits is not so strong as in the baker map case. 
However, scar functions are more robust than the other 
states, as in the previous model. 

These results show that this behavior is independent of 
the kind of noise considered, when compared to what has 
been shown in [12]. There, a different open system was 
studied (a dissipative baker map), in which the noise was 
non-generic because it acts along a preferential direction 
in phase space, corresponding to the stable manifolds of 
the classical map. Also, there is no dependence on the 
kind of map. In fact, in a more generic system in terms 
of scarring like the cat map, the same has been found. 

In order to propose an explanation for the behavior of 



scar functions, in Figs 5 and 6 we show the Wigner distri- 
butions corresponding to the first 4 steps of the evolution 
of the states used in Figs. 3 and 4, respectively. In the up- 
per panels we see the eigenstates of the maps, which have 
a complicated background structure along with localiza- 
tion on the corresponding orbit. The scar functions and 
POMs have a much simpler shape, with details that in 
general live longer than those present in the eigenstates. 
In fact, the eigenstates seem to converge to the corre- 
sponding scar functions. This is the underlying mecha- 
nism that produces a faster loss of fidelity and purity in 
the eigenstates with respect to the classically motivated 
quantum distributions. Finally, a brief discussion about 
the similar behavior of the decay of the purity and the fi- 
delity. The first quantity measures the rate of coherence 
loss, that can be seen very clearly through the Wigner 
distributions. The second one measures the correlation 
between the initial state and the evolved ones. The in- 
teraction with the environment essentially destroys the 
interference fringes whilst the dynamics distorts the ini- 
tial distributions. But all three initial states are very 
localized on a periodic orbit and its manifolds. Then the 
main initial contribution to the loss of both, purity and 
fidelity comes from the destruction of coherences. 
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Purity 

Baker map - Orbit 01 - 8 = 0. 1 




Map eigenstate (overlap = 0.828) 
Scar function 
POM 



4 5 6 

Time (t) 



Fidelity 

Baker map - Orbit 01 - 8 = 0.1 




Map eigenstate (overlap = 0.828) 

Scar function 

POM 



Figure 3: Purity and fidelity behavior for the noisy baker map. In the upper panels we show the Husimi 
representation of two initial states, localized near the period-2 orbit 01. Black corresponds to maximum probability, 
and white to minimum. Upper left panel: map eigenstate. Upper right panel: scar function. Middle panel: purity 
evolution (logarithmic scale). Lower panel: fidelity evolution (logarithmic scale). In the last two panels, black dotted 
lines correspond to the map eigenstate, green dashed lines to the POM, and red solid lines to the scar function. We 

have taken e = 0.1, N = 100. 
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Purity 



Cat map - Period 3 orbit - 8 = 0. 1 




Fidelity 



Cat map - Period 3 orbit - 8 = 0. 1 



Map eigenstate (overlap = 0.692) 
Scar function 
POM 




Figure 4: Purity and fidelity behavior for the noisy cat map. In the upper panels we show the Husimi representation 
of two initial states, localized near a period-3 orbit. Colors, patterns, scales and parameters as in Fig. 3. 
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Figure 5: Wigner distributions corresponding to the first 4 steps of evolution of a baker map eigenstate (upper 
panels), the corresponding scar function (middle panels), and POM (lower panels). We have taken the same 

parameters as those chosen for Fig. 3. 



V. SUMMARY AND CONCLUSIONS 

We have studied two complementary decoherence mea- 
sures purity and fidelity for a generic diffusive noise in 
two different chaotic maps. 

For both quantities, we have found classical structures 
in quantum mechanics that are specially stable when sub- 
jected to environmental perturbations. They are the scar 
functions, which are associated to periodic orbits and the 
stable and unstable manifolds in their vicinity. 

We have seen that quantum states constructed on clas- 
sical invariants, periodic orbits and their stable and un- 
stable manifolds, are more stable against an external 
noise than the eigenstates of the closed quantum system. 
This turns them into the most robust significant quantum 
distributions in generic dissipative maps. We conjecture 
that the scar functions will be the most robust structures 
also in general chaotic systems. 

This result has already been announced for a particu- 
lar dissipative noise in the Baker map. But in that case 
the noise was along the stable direction of the hyperbolic 



structure of the original map. Here, we have confirmed 
the same result for a general diffusive model for two dif- 
ferent maps. 

We can then say that the external noise destroys the 
stability of the quantum invariants faster than the sta- 
bility of the classical ones. This is a consequence of the 
effect of the noise as it quickly destroys quantum interfer- 
ences whilst only spreads the classical structures in phase 
space. Despite we have only shown results for e = 0.1, 
we have verified that they represent the generic behav- 
ior for a wide range of couplings. For strong couplings 
(e > 1) the noise destroys the distributions very fast, and 
for weak ones (e < 0.001) the system behaves in a similar 
way to the closed one. 

Then, we have shown that for generic maps scar func- 
tions represent the stable classical skeleton of the map 
eigenstates against environmental perturbations. 

We are currently developing the theory to quantita- 
tively explain this behavior. 
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Figure 6: Wigner distributions corresponding to the first 4 steps of evolution of a cat map eigenstate (upper 
panels), the corresponding scar function (middle panels), and POM (lower panels). We have taken the same 

parameters as those chosen for Fig. 4. 
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